ON  THE  SUPPLY  AND  DEMAND 
FOR  ENERGY 


Farihorz  Golshani  Javadi 


DUDLEY  KNOX  Llbh    . 

NAVAL  POSTGRADUATE  S_,,uul 

MONTEBEY,  CALIFORNIA  93940 


AVAL  POSTGRADUATE  SCHO 

Monterey,  California 


T 


ON  THE 

SUPPLY  AND 
FOR  ENERGY 

by 

DEMAND 

Fariborz  Golshani 

Javadi 

June  1974 

Thesis 

Ad 

visor : 

Davi 

d 

Wh 

ipple 

Approved  for  public  release;  distribution  unlimited. 


T 161523 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Date  1 

Intered) 

REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

1.     REPORT   NUMBER 

2.  GOVT  ACCESSION  NO. 

3.      RECIPIENT'S  CATALOG  NUMBER 

4.     TITLE   find  Subtitle) 

ON  THE    SUPPLY   AND   DEMAND   FOR   ENERGY 

5.     TYPE  OF   REPORT  «  PERIOD  COVERED 

Master's   Thesis 
June    1974 

6.     PERFORMING  ORG.    REPORT  NUMBER 

7.     AUTHOR(»J 

Fariborz  Golshani   Javadi 

8.     CONTRACT  OR  GRANT  NUMSERft) 

9.     PERFORMING  ORGANIZATION   NAME   ANO   AODRESS 

Naval   Postgraduate   School 
Monterey,    California      93940 

10.     PROGRAM   ELEMENT.  PROJECT,   TASK 
AREA  *   WORK   UNIT   NUMBERS 

II.     CONTROLLING  OFFICE  NAME   AND   ADDRESS 

Naval   Postgraduate   School 
Monterey,    California     93940 

12.     REPORT   DATE 

June   1974 

13.     NUMBER  OF   PAGES 

83 

14.     MONITORING  AGENCY  NAME  &   ADDRESSf"  dlllerant  Irom  Controlling  Olllce) 

Naval   Postgraduate   School 
Monterey,    California      93940 

15.     SECURITY  CLASS,  (ol  thle  report) 

Unclassified 

15«.     DECLASSIFI  CATION/ DOWN  GRADING 
SCHEDULE 

16-     DISTRIBUTION   STATEMENT  (ol  thle  Report) 

Approved   for  public   release;    distribution  unlimited. 

17.     DISTRIBUTION  STATEMENT  (ol  the  cbetract  entered  In  Block  20,  II  dltlere.nl  (rem  Report; 

18.     SUPPLEMENTARY  NOTES 

19.     KEY  WORDS  (Continue  on  reveree  tide  II  neceeemry  end  Identity  by  block  number) 

Conjunctive  management   of   energy   resources;    Analytical   Model 
for  management    of    energy    resources;    Supply    and  Demand    for   Energy; 
Empirical   and  Analytical   Study   of   Supply   and  Demand   for  Energy; 
Energy   Management . 

20.     ABSTRACT  (Continue  on  reverie  elde  II  neceeeary  end  Identity  by  fcJock  number; 

This   document    begins   with   a   review  of    the   world   energy    con- 
sumption,   it    is   shown    that    there    is    a   strong   relationship   be- 
tween   economic   growth   and   energy   consumption.       This    serves   as 
background    for    the   development    of    a   comprehensive    analytical 
model    capable   of    quantitively   evaluating   the    impact    of    energy 
related   decisions.      The   general   model   developed    is   descriptive 
of    exploration,    extraction,    storage,    import    and  processing   of 

DD  ,^NRM7J  1473 
(Page    1) 


EOlTlON   OF    1   NOV  65  IS  OBSOLETE 
S/N    0102-014-  6601   I 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wnen  Dele  entered) 


CtCUWITY   CLASSIFICATION   OF    THIS  P*GEf*)im  D»l«  Enl»rvJ' 


Block    20    -    ABSTRACT    (Cont.) 

energy   resources.      The   analytical  model   also   takes   into   account 
the   relationship  between   these   aspects   of   energy  production 
and   the  market   prices   of   energy   resources. 


DD      Form       1473  (BACK) 

.  1  Jan  73 

S/N    0102-014-GG01  security  classification  of  this  pageo^i»"  D«m  em»'»<i> 


ON  SUPPLY  AND  DEMAND  FOR  ENERGY 


by 


Farihorz  ,Golshani  Javadi 
Lieutenant,  Iranian  Navy 
B.S.,  Italian  Naval  Academy,  1965 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  MANAGEMENT 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
June  1974 


ft' 


ABSTRACT 

This  document  begins  with  a  review  of  the  world  energy 
consumption,  it  is  shown  that  there  is  a  strong  relationship 
between  economic  growth  and  energy  consumption.   This  serves 
as  background  for  the  development  of  a  comprehensive  analyt- 
ical model  capable  of  quantitively  evaluating  the  impact  of 
energy  related  decisions.   The  general  model  developed  is 
descriptive  of  exploration,  extraction,  storage,  import  and 
processing  of  energy  resources.   The  analytical  model  also 
takes  into  account  the  relationship  between  these  aspects  of 
energy  production  and  the  market  prices  of  energy  resources. 
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I.   INTRODUCTION 

One  of  the  most  important  issues  in  the  world  today  is 
the  soaring  prices  of  energy  resources,  this  has  caused  con- 
siderable anxiety  in  the  energy  consuming  nations  concerning 
their  ability  to  maintain  economic  growth  and  financial 
stability. 

For  example,  in  his  second  economic  report,  President 

Nixon  cautioned  that: 

"The  new  problems  we  face  are  of  such  enormity 
that  there  may  be  a  temptation  to  delay  further 
progress  toward  trade  and  monetary  reform.   The 
drastic  increase  in  oil  prices  will  have  a  signif- 
icant impact  on  both  the  domestic  economies  of  all 
nations  and  international  economic  relationships . [1] 

According  to  United  Nations  Statistical  Office,  the  total 
energy  consumption  in  the  world  increased  by  4%  in  1971. 
The  world  consumption  of  electricity  increased  by  6%  over 
previous  year;  world  production  of  crude  petroleum  increased 
by  5.5%;  the  world  exports  of  natural  gas  increased  by  36% 
over  1970.   The  nuclear  power  accounted  for  2%  of  total  en- 
ergy produced,  which  was  an  increase  of  38%  over  1970.   On 
the  other  side  the  use  of  solid  fuels  declined  by  9%  of 
total  world  energy  consumption.  [2] 

Government  revenues  in  oil  producing  nations  have  in- 
creased considerably  due  to  the  recent  increases  in  the  price 
of  crude  oil.   In  fact  the  posted  price  of  Persian  Gulf  light 
crude  oil  was  increased  from  U.S.  $1.80  per  barrel  on  June  1, 
1973  to  U.S.  $11.50  per  barrel  on  January  1,  1974  by  OPEC 
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(organization  of  Petroleum  Exporting  Countries),  consequently 

producing  nations  are  accumulating  more  reserves  and  their 

surplus  could  rise  from  $10  billion  in  1973  to  $60  billion 

in  1974.  [3] 

On  the  other  side,  consuming  nations,  due  to  the  higher 

oil  prices  would  have  to  have  current  account  deficits.   In 

their  conference  in  Washington  the  consuming  nations  agreed 

that:    "...the  present  situation,  if  continued  could 
lead  to  a  serious  deterioration  in  income  and 
employment,  intensify  inflationary  pressures, 
and  endanger  the  welfare  of  nations.   And  a 
serious  setback  to  the  prospect  for  economic 
development  of  developing  countries  will  re- 
sult from  their  additional  energy  costs  at 
current  oil  prices."  [4] 

The  production  and  consumption  of  energy  has  significant  im- 
pact on  the  economy  and  financial  position  of  a  country. 
Hence  in  the  last  few  years,  many  studies  have  been  conducted 
to  evaluate  the  impacts  of  policy  options  such  as  price  reg- 
ulations, import  regulations,  environmental  controls,  etc. 
and  several  of  these  studies  have  used  advanced  analytical 
techniques. 

In  order  to  be  able  to  trace  the  impact  of  policy  changes 
and  to  be  able  to  develop  an  analytical  model  for  the  con- 
junctive management  of  the  energy  sector  of  economy,  it  is 
essential  to  have  an  understanding  of  some  of  the  complex 
interrelationships  in  the  broad  energy,  demand  and  pricing 
situation . 

This  is  the  primary  objective  of  this  study.   By  conduct- 
ing a  conprehensive  analysis  of  the  world's  energy  consump- 
tion, energy-GNP  relationship,  and  the  existing  literature, 


an  analytical  model  has  been  developed.   The  usefulness  of 
this  study  for  future  efforts,  in  analyzing  the  consequences 
of  alternative  energy  policy  options,  is  covered  in  the  con- 
clusion. 
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II.   AN  EMPIRICAL  STUDY  OF  THE  ENERGY-GNP  RELATIONSHIP 

A.   INTRODUCTION 

"Economic  history  attests  to  the  crucial  role 
played  by  the  consumption  of  energy  in  advancing 
the  material  well-being  of  mankind.   The  industrial 
revolution  and  the  growth  of  industry  in  the  nine- 
teenth century  are  almost  synonymous  with  the  sig- 
nificant contribution  of  coal  to  the  development 
of  the  iron  and  steel  industry,  to  railways,  and 
to  factory  mechanization.   In  the  twentieth  century, 
electrification  and  motorized  transport  served  to 
modify  the  process  and  many  ways  helped  economic 
progress.   'High-energy  civilization1  has  by  now 
become  a  commonplace  characterization  of  economi- 
cally-advanced 20th  Century  society.   The  growth 
of  population  and  the  levels  of  prosperity  achieved 
have  depended  upon,  and  resulted  in  (among  other 
critical  factors)  the  consumption  of  energy  in 
large  amounts.   To  the  extent  it  materializes,  con- 
tinued demographic  and  economic  growth  will  un- 
doubtedly involve  the  steadily  rising  use  of  ener- 
gy." [5] 

It  therefore  becomes  important  to  give  particular  attention 
to  this  relationship.   At  a  minimum,  the  topic  calls  for  an 
interlinking  of  physical  energy  data  and  of  the  nations'  in- 
comes and  products  (GNP)  accounts. 

Such  a  study  could  yield  a  number  of  important  results: 
A  better  understanding  of  the  behavior  of  energy 
consumption  relative  to  gross  national  product. 
It  would  attempt  to  answer  the  question:  On  the  aver- 
age what  levels  of  energy  consumption  are  required 
for  given  levels  of  GNP.  [6] 
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B.   DEFINITION  OF  THE  MODELS 

We  now  turn  to  the  specification  of  a  basic  model  relat- 
ing GNP  and  energy  consumption. 

•  In  the  first  model  the  assumption  is  made,  that  a  linear 
relationship  exists  between  per  capita  energy  consumption 
and  GNP  per  capita,  then: 

(1)  E  =  a+8G+e 

where;   E  =  energy  per  capita,  (Dependent  variable) 
G  =  GNP  per  capita,  (Independent  variable) 
a  =  constant 
3  =  the  coefficient  representing  the  propensity  to 

consume  energy,  (in  physical  units) 
e   =  stochastic  disturbance. 
Next  the  assumption  is  made,  that  an  exponential  rela- 
tionship or  logarithmic  relation  exists  between  these  vari- 
ables, then: 

(2)  E  =  yG6I  ,  or 

(3)  In  E  =  lny  +  61nG  +  lne 
where;  y   =   lny  =  constant. 

6  =  coefficient  representing  the  elasticity  of  energy 
consumption  with  respect  to  GNP 

e  =  lne  =  stochastic  disturbance 
The  models  represented  by  equations  (1)  and  (3)  were  applied 
for  a  number  of  samples  of  data  corresponding  to  OECD  Coun- 
tries, (Organization  for  Economic  Cooperation  and  Develop- 
ment).  The  results  of  the  regression  analysis  are  presented 
in  section  D,  and  the  summary  of  principal  statistical  re- 
sults appear  in  Table  I. 
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C.   DESCRIPTION  OF  DATA 

The  data  used  here  are  collected  from  United  Nations 
Statistical  yearbooks  and  Statistics  of  Energy  for  member 
countries  of  OECD  listed  in  Appendix  C.   The  selection  of 
data  was  made  because  of  the  special  nature  of  the  problem 
which  requires  detailed  statistics  of  GNP  in  constant  prices, 
with  corresponding  values  for  population  and  energy  consump- 
tion.  The  only  available  data  with  the  appropriate  charac- 
teristics were  those  of  OECD  countries  and  only  for  selected 
years,  the  specific  data  which  is  used  is  described  in  fol- 
lowing lines. 

1.  Population 

Data  are  taken  from  OECD  Main  Economic  Indicators, 
and  is  measured  in  million,  for  years  1960,  1965,  1966,  1967, 
1968,  1969,  1970  and  1971,  total  of  174  observations  (22 
countries  multiplied  by  8  years  minus  two  observations  for 
years  1966,  1967  of  Australia,  not  being  used,  because  of 
lack  of  corresponding  values  for  energy  consumption  per  capi- 
ta during  these  two  years.) 

2 .  Energy  Consumption  per  Capita 

Data  are  taken  from  United  Nations  Statistical  year- 
books; 1960,  1965,  1966,  1967,  1968,  1969,  1970  and  1971. 
It  is  measured  in  kilograms  coal  equivalent  per  capita. 
There  are  174  observations. 

3 .  Gross  National  Product 

Data  are  taken  from  OECD  Main  Economic  Indicators. 
GNP  is  measured  in  million  U.S.  dollars,  at  1963  prices  and 
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1963  exchange  rates.   The  GNP  figures  are  next  divided  by 
corresponding  values  of  population  to  determine  the  GNP/ 
capita,  for  the  years  1960,  1966,  1967,  1968,  1969,  1970  and 
1971  with  a  total  of  174  observations. 
4.   By  Sources  of  Energy 

All  data  are  taken  from  "Statistics  of  OECD  countries 
for  Energy,  1954-1968,  Paris  1970"  for  the  years  1966,  1967 
and  1968,  and  divided  by  corresponding  population  figures 
(Appendix  C). 

a.  Crude  Petroleum 

This  is  measured  in  thousands  metric  tons  with  a 
total  of  60  observations  (20  countries  for  the  three  years, 
1966,  1967  and  1968). 

b.  Motor  Gasoline 

This  is  measured  in  million  metric  tons  with  a 
total  of  63  observations  (21  countries  for  the  three  years, 
1966,  1967  and  1968). 

c.  Aviation  Fuel 

This  is  determined  by  sum  of  the  aviation  gaso- 
line and  jet  fuel.   It  is  measured  in  million  metric  tons 
with  a  total  of  63  observations  (21  countries  for  the  three 
years,  1966,  1967  and  1968). 

d.  Hard  Coal 

This  is  measured  in  million  metric  tons  with  a 
total  of  63  observations  (21  countries  for  the  three  years, 
1966,  1967  and  1968). 
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e.  Kerosine 

This  is  measured  in  million  metric  tons  with  a 
total  of  63  observations  (21  countries  for  the  three  years, 
1966,  1967  and  1968). 

f.  Gas  Oil 

This  is  measured  in  million  metric  tons  with  a 
total  of  63  observations  (21  countries  for  the  three  years, 
1966,  1967  and  1968). 

g.  Natural  Gas 

This  is  measured  in  billion  cubic  meters  at  4200 
KCal.,  with  a  total  of  33  observations  (11  countries  for  the 
three  years,  1966,  1967  and  1968). 

D.   STATISTICAL  RESULTS 

For  the  application  of  the  models  a  general-purpose  sta- 
tistical package,  based  on  least-square  method  (SNAP/IEDA) 
has  been  used.   The  results  of  the  analysis  are  as  follows 
(with  standard  errors  in  parenthesis), 
where;   R2  =  Square  of  correlation  coefficient. 
R2  =  Adjusted  R2. 

F   =  The  ratio  of  the  variance  of  the  residuals  of 
the  dependent  variable  before  regression  and 
the  variance  of  the  residuals  of  that  variable 
after  the  regression. 
NOBS  =  Number  of  observations. 
M   =  The  mean  of  variables. 
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K  =  Elasticity  of  consumption.   When  the  linear  mod- 

i  •      a   v         9E   MG         3E  .   ... 
el  is  used  K  =  -gg  rrp,  where  -s-=  is  the  regres- 
sion coefficient,  MG  and  ME  are  the  means  of 
the  independent  variable  and  dependent  variable 
respectively.   When  the  constant  elasticity 
model  is  used,  K  equals  the  coefficient  5. 
Energy/Capita  and  GNP/Capita 

a.  Linear  Model  (Fig.  1) 

E  =  -584.905  +  2.625G 

(.093) 

R2  =  .823,   F  -  797.1,   NOBS  =  174, 

ME  =  3722.8621,   MG  =  1641.6658, 

k  -  I"!  -   2-625    «  iflHfi  -   2.625    .    .441 

K   =    1.158,      I2    =    .822 

b.  Constant  Elasticity  Model  (Fig.  2) 

InE  =  -0.902  +  1.225  InG 

(.031) 

R2  -  .903,   F  =  1603.4,   NOBS  -  174, 

K  =  1.225,   R2  =  .902 

Crude  Petroleum/Capita  and  GNP/Capita 

a.   Linear  Model  (Fig.  3) 

Ep  =  269.059  +  0.709G 

(0.269) 

R2  =  .107,   F  =  6.9,   NOBS  =  60, 

MEp  =  1401.8915,   MG  =  1597.0630, 

1597.0630 
K  -  .709  x  -— __— 

K  -  .709  x  1.14  =  .81,   R2  -  .092 
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b.   Constant  Elasticity  Model  (Fig.  4) 

InEp  =  -0.334  +  0.998  InG 

(.126) 

R2  =  .52,   F  =   62.8,   NOBS  =  60 

K  =  .998,   R2  =  .488 

Motor  Gasoline/Capita  and  GNP/Capita 

a.  Linear  Model  (Fig.  5) 

Eg  =  -172.334  +  0.258G 

(.021) 

R2  =  .705,   F  =  145.9,   NOBS  =  63, 

MEg  =  244.5014,   MG  =  1616.1303, 

K  "  ,258    244.5014 

K  =  .258  x  6.610  =  1.705,   R2  =  .700 

b.  Constant  Elasticity  Model  (Fig.  6) 

InEg  -  -4.421  +  1.321  InG 

(.061) 
R2  =  .885,   F  =  470.6,   NOBS  =  63, 

K  =  1.321,   R2  =  .883 

Aviation  Fuel/Capita  and  GNP/Capita 

a.   Linear  Model  (Fig.  7) 

Ea  =  -23.644  +  0.045G 

( .006) 

R2  =  .451,   F  =  50.1,   NOBS  =  63, 

MEa  =  48.4012,   MG  =  1616.1303, 

K  =   045  x  1616.1303 
•U^°    48.4012 

K  =  .045  x    33.390  =  1.503,   R2  =  .442 
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b.   Constant  Elasticity  Model  (Fig.  8) 

InEa  =  -4.304  +  1.061  InG 

(.166) 

R2  =  .402,   F  =  40.9,   NOBS  =  63, 

K  =  1.061,   R2  =  .392 

5.  Hard  Coal/Capita  and  GNP/Capita 

a.  Linear  Model  (Fig.  9) 

Eh  =  69.885  +  0.470G 
(.123) 
R2  =  .192,   F  =  14.5,   NOBS  =  63, 

MEh  =  828.7752,   MG  =  1616.1303, 

K  =   470  x  1616.1303 
K    '47°  X  838.7752 

K  =  .470  x  1.95  =  .92,   R2  =  .179 

b.  Constant  Elasticity  Model  (Fig.  10) 

InEh  =  -.917  +  .957G 
(.267) 

R2  =  .174,   F  =  12.9,   NOBS  =  63, 

K  =  .957,   R2=.160 

6.  Kerosine/Capita  and  GNP/Capita 

a.  Linear  Model  (Fig.  11) 

Ek  =  6.973  +  .014G 
( .005) 
R2  =  .121,   F  =  8.4,   NOBS  =  63, 

MEk  =  29.6929,   MG  =  1616.1303, 

K  -   014  x  161S-1303 
K  ~  -014  X  29.6939 

K  =  .014  x  54.43  =  .762;   R2  =  .107 

b.  Constant  Elasticity  Model  (Fig.  12) 
R2  =  .06,   R2  =  .045 
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7.  Gas  Oil/Capita  and  GNP/Capita 

a.  Linear  Model  (Fig.  13) 

Eo  =  229.796  +  .452G 
(.065) 

R2  =  .443,   F  =  48.6,   NOBS  -  63, 

MEo  =  960.6031,   MG  =  1616.1303 

1616.1303 
K    "4DZ    960.6031 

K  =  .452  x  1.682  =  .760,   R2  =  .434 

b.  Constant  Elasticity  Model  (Fig.  14) 

InEo  -  -1.462  +  1.121  InG 

(.071) 

R2  =  .805,   F  =  251.7,   NOBS  =  63, 

K  =  1.121,   R2  =  .802 

8.  Natural  Gas/Capita  and  GNP/Capita 

a.  Linear  Model  (Fig.  15) 

En  =  -2013.059  +  1 . 708G 

(.184) 

R2  =  .736,   F  -  86.4,   NOBS  =  33, 

MEn  =  992.0987,   MG  =  1759.8145, 

1759.8145 
K  -  1.708  x  992  0987 

K  =  1.708  x  1.77  =  3.03,   R2  =  .727 

b.  Constant  Elasticity  Model  (Fig.  16) 

InEn  =  -27.099  +  4.392  InG 

(.609) 

R2  -  .626,   F  =  52.0,   NOBS  =  33, 

K  =  4.392,   R2  =  .614 

The  results  of  the  empirical  analysis  are  summarized  in 

Table  I. 
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E.   INTERPRETATION  OF  THE  RESULTS 

By  referring  to  the  summary  of  results  presented  by  Table 
I,  we  can  reach  the  following  conclusions: 

1)  The  logarithmic  model  was  a  better  fit  to  the  data 
on  energy  per  capita,  crude  petroleum,  gasoline,  and  gas  oil. 
While  the  simple  linear  model  was  a  better  fit  to  data  on 
aviation  fuel,  hard  coal,  kerosine  and  natural  gas. 

2)  There  is  a  higher  correlation  between  GNP  and  aggre- 
gate energy  consumption  than  between  GNP  and  different  ener- 
gy sources.   This  could  be  explained  by  the  fact  that  each 
member  country  of  the  OECD  consumes  different  sources  of 
energy  as  the  primary  mixture  for  energy  production.   There 
are  a  number  of  economic  and  technological  factors  affecting 
the  mix  of  the  primary  energy  sources  within  a  country.   It 
is  therefore  necessary  to  pay  particular  attention  to  the 
special  factors  underlying  the  historic  interfuel  composition 
of  total  energy  consumption  and  to  the  policy  tools  avail- 
able to  alter  the  composition. 
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III.   THE  ANALYTICAL  STUDY  OF  ENERGY  PRODUCTION  AND  MANAGEMENT 

A.   INTRODUCTION 

There  are  considerable  efforts  taking  place  in  the  private 
and  public  sectors  of  the  economics  profession  in  the  recent 
years,  to  satisfy  the  need  for  a  quantitative  economic  theory 
of  energy  conservation  and  management.   Such  efforts  provide 
a  point  of  departure  for  the  development  of  analytical  mod- 
els which  may  be  usefule  in  projecting  future  energy  require- 
ments and  the  short-term,  and  the  long-term  adequacy  of 
supplies  to  meet  them. 

The  energy  situation  comprises  a  wide  range  of  interre- 
lated socio-economic  and  institutional  factors.   These  in- 
clude: 

-  Balance  of  payments. 

-  National  security  and  stockpiling. 

-  Environmental  costs  and  user  costs. 
Protection  of  energy  related  industries. 

-  Resource  production  (exploration,  development,  ex- 
traction, processing,  refining)  and  relative  facil- 
ities. 

Transportation  (means,  modes  and  facilities). 

-  Marketing  and  Distribution  systems. 
Substitutabilities  among  different  energy  sources. 

-  Development  of  new  sources  of  energy. 


22 


Problem  involving  risks  and  uncertainties  in  terms 
of  domestic  and  foreign  investments  in  exploration 
as  well  as  for  import  supplies,  prices,  and  uncer- 
tainties related  to  technologies. 
-  Interdependencies  among  the  stages  of  activity  and 
the  policy  alternatives. 

Conjunctive  management  at  national,  regional,  and 
state  level  versus  an  individual  firm  or  a  single 
industry. 
The  purpose  of  this  chapter  is  to  present  a  summary  of 

the  literature  concerning  the  production  and  management  of 

the  natural  resources. 

B.   RESULTS  OF  THE  ECONOMIC  THEORY  OF  EXHAUSTIBLE  RESOURCES 
The  theory  of  exhaustion  considers  the  losses  firms  may 
incure  if  they  do  not  consider  the  resource  availability  in 
future  generations. 

Results  indicate  that  different  structure  of  industry  may 
produce  different  output  patterns,  for  example  private  owner- 
ship may  lead  to  decisions  different  from  those  of  central- 
ized management.   In  addition,  the  time-path  of  exploitation, 
when  the  costs  are  affected  by  cumulative  production  is  dif- 
ferent from  conditions  where  the  increasing  costs  due  to 
cumulative  output  are  ignored. 

The  general  format  of  the  models  presented  under  Exhaus- 
tion theory  consist  of  maximization  of  present  value  of 
future  profits  subject  to  limited  supply  of  resources.   In 
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following  section  simple  examples  of  these  models  have  been 
presented. 

C.   SIMPLE  CASES  IN  THEORY  OF  EXHAUSTION 
1.   Case  (1) 

Considering  the  case  where  a  private  distributor  of 
fuel  is  hoarding  a  fixed  quantity  of  it  for  a  variable  length 
of  time  without  occurring  additional  costs,  and  the  selling 
price  of  fuel  increases  by  a  known  function  P  =  p(t);  as 
time  goes  on.   What  would  be  the  optimum  time  for  him  to  sell 
stored  supply  knowing  that  he  can  sell  all  quantity  Q  avail- 
able? [7] 

It  is  assumed  that  interest  is  added  continuously  at 
a  constant  rate  of  r%. 

If  (y)  denotes  the  present  value  of  the  fuel  sold 
after  (t)  years,  then 

(1)  y  =  Q-P(t)-e~rt 

The  condition  for  optimum  time  are  given  by: 

dt    U  '    dt*   ° 
Knowing  that  the  log  y  is  a  "monotonic  transformation"  of  y. 
(This  means  that  whenever  y  is  increasing  log  y  is  also  in- 
creasing, and  whenever  y  is  falling  its  logarithm  is  also 
falling.)   Therefore,  the  point  corresponding  to  the  maximum 
of  y  is  also  maximum  of  log  y.   Then 

(2)  log  y  -  log  P(t)  +  log  Q  -  rt,   or 

(3)  log  y  =  F(t)  +  log  Q  -  rt 
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where    F(t)  =  log  P(t).   Then: 

fEClo«,)-i$-f(t)-r 


d 


zdog  y)  -  Uk  S>  -  I  ^  -  ir(^)2  ■  At) 


dfv   &  J  '    dtvy  dt'    y  dt z   y2Vdt 

The  condition  for  optimum  selling  time  are  thus 

F"(t)  <  0   and   F  (t)  =  r 

From  Eq.     (3)    the    first    condition    is    that 

i 
P'(t)   -  ftdog  P(t))    =  ^T)    =   r 

Where  gM  determines  the  rate  of  growth  of  selling  price 
at  period  (t),  is  equal  to  the  constant  rate  of  market  inter- 
est r.   By  the  second  condition 

2 


"       d  d   Pft^ 

F  (t)  -JfTUog  P(t))  =^ltl)}  <  ° 


which  implies  the  rate  of  growth  for  selling  price  must  be 
decreasing  at  the  optimum  time. 

Result:   If  the  rate  of  increase  for  selling  price 
is  constant  and  equal  to  the  market  rate  of  interest  over  the 

time;  then  owner  would  be  indifferent  whether  he  sells  his 

rt 
reserves  at  price  Pn  now  or  at  a  price  P=Pne   after  time  t. 

This  could  be  shown  from  Eq .  (1).   The  present  value  would 

be:      y  =  QPe"rt 

rt 
where  P=p0e    implies 

(4)      y  =  QP0ert.e"rt  =  QPQ 


Appendix  A,  Section  1. 
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2.   Case  (2) 

Now  considering  the  case  under  free  competition  when 

the  rate  of  increase  for  selling  price  is  equal  to  the  mar- 

rt 
ket  rate  of  interest.   The  formula  P=Pne    fixes  the  relative 

prices  at  different  times  and  the  owner  would  be  indifferent 

about  a  unit  of  fuel  being  sold  now  or  in  the  future.   (P 

being  interpreted  as  the  net  price  received  after  paying  the 

costs  of  extraction  and  placing  upon  the  market.)  [8]   And  Q 

being  the  total  amount  of  supply  available,  and  the  quantity 

q  distributed  at  any  given  time  is  a  continuous  function  of 

price  and  time,  if  t=T  being  the  time  of  final  exhaustion, 

it  could  be  written: 

(5)  q  =  q(P,t)    when  t=0,  P=PQ;  then  q=qQ- 

rT 
and  when  t=T,  P=P  e   and  q=0.   Then 

T  T 

(6)  Q  =  /  q(P,t)dt  =  /  q  (P  ert,t)dt 

0  0   u   u 

The  question  will  be  the  determination  of  t=T  the 
time  of  final  exhaustion  to  maximize  the  present  value  of 
his  revenue.   Since  at  t=T,  q=0  it  implies: 

(7)  q  =  q(P.t)  =  q(P0erT,T)  =  0 

and  consequently 

T  T 

y  =  /  q(P,t)-P-e~rtdt  =  /  q(P , t ) -Pert •e~rtdt 
0  0 


=  /  q(P,t)-P  dt  =  P   /  q(p,t)dt 
0  U       U  0 

From  Eq.  (6)  follows  y=P0Q  where  P_  and  Q  are  both  constant 
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Result:   In  the  case  of  free  competition,  when  also 
the  rate  of  price  increase  is  constant  and  equal  to  market 
rate  of  interest,  the  owner  is  indifferent  to  q  quantity 
supplied  at  any  time  t  from  Q  his  total  reserve  and  the  T 
final  time  of  exhaustion. 

D.   A  GENERAL  CASE  IN  THEORY  OF  EXHAUSTION 

Considering  the  case  under  monopoly  when  a  certain  com- 
modity is  available  in  fixed  supply  Q  with  the  rate  of  pro- 
duction q=q(t)  over  a  time  period  from  t0  to  t„  [9,  10]. 
Production  then  is  subject  to  the  constraint. 

t2 

(8)  /  q(t)dt  <  Q 

to 
The  firm  would  plan  to  maximize  the  present  value  of  its 
profits 

t2 

(9)  y  =  /  n[q(t),t]-e  rtdt 

to 

where  y  is  net  present  value;  II  total  profits;  q(t)  output; 
r  the  continuous  interest  rate. 

The  question  involves  the  selection  of  the  optimal  times 
for  starting  and  completing  output  and  the  optimal  output 
pattern.   The  solution  implies  resolving  the  Lagrangean  Equa- 
tion. 

t2  -rt       t2 

(10)  h  =    f   n[q(t),t]-e   udt  -  X    f   q(t)dt-Q  . 

to  to 

where  A  is  the  Lagrangean  multiplier;  and  Q  total  availabil- 
ity. 
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The  optimal  production  pattern  is  determined  by  the 

2 

Euler  equation  of  the  Calculus  of  Variation,   and  the  bound- 
ary conditions  are  the  attaining  stationary  value  at  t 0  and 
t2  (Eq.  10).   The  Euler  equation  requires  that 

oLi    9  f  9L  ^  _  « 
"3q  "  3tllq; 

And  because  h~~^   does  not  appear  in  Eq.  (10),  it  therefore 

does  not  effect  profit,  so  the  Euler  equation  reduces  to 

/1-.N     3L    -rt  r29II  9q..    ,     _ 

(11)  "Sq  =  e    '  To"  ^tdt  -  Xq  =   ° 

t  o 

By  substitution  of  dq  for  -r^fdt ;  Eq.  (10)  could  be  written  as 

(12)  ||  ,  e-rt  /2Mdq  -  xc,   or 

t  0 

2  911        rt 

(13)  /  -k— dq  =  Xqe  or 
.  dq 

t  o 

(14)  -^  =  Xe 

By  introducing  the  concept  of  marginal  profit,  the  revenue 

(R)  minus  Cost  (C)  would  determine  the  profit.   Using  letters 

^R 
P  for  price;  AC,  average  costs;  <J>,  average  profits;  MR  =  -r— , 

dq 

9C  9n 

marginal  revenue;  MC  =  -r— ,  marginal  cost;  and  M   =  — ,  mar- 

dq  dq 

ginal  profit.   Then  Eq .  (14)  can  be  written  as 

(15)  MII(t)  =  -|^  =  MR  -  MC  =  Aert 

9q 

9 2  n 

The  conventional  second  order  equations  - — j-  <  0  and 

92R    92C   ,       . 
j^t   <  -^r   also  apply. 


2 

Appendix  B, 
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In  fact  n  =  R  -  C  =  P(q)-q  -  C(q), 

Mn=l^  =  l|"lf  =  P(q)  +  (1*P(C1)  ~  C^>  and' 

(16)  1^  =  f^#  -  f£   =  P(q)  +  P(q)  +  q-P(q)  -  C(q) 

where  in  general  P(q)  being  the  demand  curve  implies  its 

•  it 

slope  P(q)  <  0  and  P(q)  _>  0.   Considering  increasing  cost 

it       ii 
industry,  it  then  implies  C(q)  >  P(q)  >  0;  consequently  mak- 
ing Eq.  (16)  negative  or 

d2Il    „  _      32R    32C  ,    n  32R  „  32C 

■z — z   <  0,  or  -r — r  -  -5 — t   <  0,  or  ■* — T   <   -5 — r    • 
3q  8q     3q  9q     dq 

Under  free  competition,  P  =  MR  so  Eq.  (15)  becomes 

(17)  P  -  MC  =  Xert 

and  with  constant  cost,  this  reduces  further  to 

(18)  P  -  AC  =  Xert 

Result:   From  Eq.  (15)  it  derives: 

(19)  n  =  Aq(t)ert 
dividing  it  by  q(t)  results: 

(20)  <f(t)  =  Xert 

since  Eq .  (15)  also  hold  at  time  t2,  the  optimum  production 
path  is  where  marginal  profit  equals  average  profit.   Since 

H(t)  =  <Kt)-q(t),  then  Mll(t)  =  <|>(t)  +  (f^)q(t). 

The  equality  MII(t)  =  4>(t)  implies  q(y^)  =  0  which  occurs  if 
q  or  tt^  is  zero;  the  firm  either  produces  nothing  at  t2  or 
produces  the  output  that  maximizes  the  average  profits  at  t2 
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E.   RESULTS  OF  ECONOMIC  MODELS  OF  PRODUCTION  IN  RESOURCE 

INDUSTRIES 

The  theory  of  exhaustion  has  been  developed  without  ex- 
plicitly considering  the  role  of  investment  associated  with 
resources,  the  interaction  among  production  rates,  invest- 
ment rates,  technological  changes  and  exploration  expendi- 
tures. 

In  recent  studies  efforts  have  taken  place  to  develop 
economic  models  to  include  these  factors.   According  to  these 
studies,  the  stages  of  activity  necessary  for  the  production 
of  energy  commodities  involve  numerous  interrelated  aspects. 
For  example,  the  optimal  extraction  rate  of  a  given  resource 
and  associated  investments  in  capital  stock  are  related  to 
exploration  and  importation  policies  of  the  resources  as 
well  as  the  substitutability  among  the  resources.  [12] 

The  supply  of  particular  resource  should  be  obtained 
from  domestic  production  and  imports  until  marginal  import 
costs  equal  the  marginal  costs  of  extracting  domestic  stocks, 
and  be  allocated  between  current  use  and  stockpiles  for  fu- 
ture consumption,  until  the  marginal  contribution  of  resource 
to  the  current  use  equals  its  marginal  contribution  to  future 
use.  [13] 

By  considering  the  user  costs  defined  by  A.  Scott  as  "the 
present  value  of  future  profits  foregone  by  a  decision  to 
produce  a  unit  of  output  today."   The  use  of  imports  in  cur- 
rent production  would  reduce  the  rates  of  domestic  extraction 
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and  consequently  the  user  cost,  or  imports  could  be  used  to 
increase  the  stockpiles  to  reduce  the  future  scarcity  of  the 
resource.  [14] 

The  above  brings  to  attention  the  costs  which  accompany 
the  import  restrictions  and/or  stock  controls. 

The  equality  of  the  user  costs  specially  among  substitut- 
able  resources  are  of  special  interest  in  the  conjunctive 
management  of  energy-related  resources. 

An  important  question  arises  in  terms  of  the  impact  of 
subsidies,  depletion  allowances  granted  to  one  resource,  on 
the  optimal  rates  of  exploration  expenditures,  importation 
and  production  of  other  resources.  [15] 

Of  particular  attention  are  the  social  attitudes  towards 
the  generation  of  pollutants,  and  the  inclusion  of  pollution 
abatement  costs  in  the  processing  costs  and  production  costs, 
which  may  be  imposed  in  terms  of  social  costs. 

Considering  the  problems  associated  with  petroleum  pro- 
duction and  reservoir  management,  the  results  suggest  that 
the  amount  of  pressure  which  exists  in  the  reservoir,  and 
the  time-path  of  pressure  depletion,  are  two  influential  fac- 
tors on  the  production  rate  and  the  amount  of  oil  recover- 
able [16] 

It  should  be  noted  that  in  many  cases  pressure  may  be 
maintained  or  increased  artificially  by  incurring  additional 
expenditures,  which  implies  that  the  recoverable  stock,  as 
well  as  production  rates,  generally  will  depend  strongly  on 
investment.   This  dependence  is  not  merely  on  a  cumulative 
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investment,  but  also  on  the  time-path  of  investment.  Thus 
the  impact  of  rapid  production  rates,  in  the  face  that  the 
recoverable  stock  is  rate  dependent,  may  be  offset  to  some 
extent  by  investments.  [17] 
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IV.   DEVELOPMENT  OF  AN  ANALYTICAL  MODEL 

A.   INTRODUCTION 

Because  of  increasing  concern  over  the  management  of 
"Energy  Sector"  of  economy  as  whole,  a  comprehensive  analyti- 
cal model  is  required,  to  be  able  to  translate  and  quantita- 
tively evaluate  the  impact  of  energy  related  procedures. 

In  the  three  preceding  chapters  a  brief  review  of  litera- 
ture, past  and  current  works  was  performed,  to  give  an  under- 
standing of  some  of  the  complex  interrelationships  in  this 
sector. 

In  specifying  the  scope  and  coverages  of  this  model,  a 
number  of  requirements,  could  be  identified,  to  be  covered 
by  the  comprehensive  model.   These  requirements  could  be 
listed  as  follows:  [18] 

The  model  should  cover  all  energy  forms. 

-  The  model  should  consider  supply  and  demand  inter- 
actions or  the  mechanism  for  balancing  supply  and 
demand. 

-  The  model  should  take  into  account  interfuel  substi- 
tution. 

The  model  must  take  into  consideration  the  impact  of 
price  and  other  variables  on  .the  demand  and  supply  of 
the  specific  fuels. 

-  The  model  must  consider  the  impact  of  technological 
changes  on  both  the  supply  and  demand  of  fuels  and 
energy. 
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The  model  must  be  disaggregated  at  the  regional  level. 

The  model  must  be  able  to  use  the  available  data  for 

future  analyses. 
Note  that  the  environmental  consequences  have  been  left  out 
in  this  study. 

The  economic  model  prepared  by  R.  Cummings  and  David 
Whipple   [19],  is  an  initial  attempt  at  a  comprehensive  model 
of  the  discovery,  extraction,  importation,  and  stockpiling  of 
primary  energy  resources,  integrating  the  production  of  final 
energy  commodities,  it  is  believed  its  further  expansion  and 
refinements,  will  make  it  suitable  for  use  in  constructing 
an  analytical  model,  presented  in  next  section,  feasible  to 
the  conjunctive  management  problems  of  energy  sector,  by 
covering  the  above  mentioned  requirements. 

B.   THE  MODEL 

By  assuming  the  existence  of  "j"  resource  fields  in  which 
any  of  "R"  resources  are  extracted  from  domestic  deposits 
and/or  stored  in  stockpiles.   For  each  field  j=l, 2 , . . . ,J;  the 
intertemporal  changes  in  stocks  of  domestic  deposits  of  the 
resource  r,  r=l,2,...,R,  are  described  by  the  following: 

(1)   Xt+1  =  X^-W1.,  X°.  given,  for  all  r=l,...,R; 
v      rj      rj   rj'   rj  '      ' 

j=l, . . . ,J;  t=l,2 T. 

By  (1),  the  sotck  of  resource  r  in  field  j  at  the  end  of 

period  t,  "X  .  "  equals  stocks  at  the  beginning  of  the  period 

"X  ."  minus  extraction  of  r  during  t,  "W  ."  which  is  a 
rj  rj 

function  of: 
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(2)  W*  .  =  <f>  .(n* .,  p*  .,  Q*), 
v     rj    Ywjv  rj'  yrj'    ^r" 

where;   n   =  The  excess  of  "P  " ;  the  market  price  of  re- 
'    rj  r  * 

wt 
source  r,  over  "C  ."  the  extraction  cost  of 

rj 

resource  r  in  field  j,  during  period  t.  [20] 

p  .  =  The  excess  of  P   over  "C  ."  the  importation 
rj  r        rj        * 

cost  of  resource  r  in  field  j ,  during  period 
t.  [21] 
Q   =  The  quantity  of  resource  r,  demanded  in  period 
t,  which  is  a  function  of  P  ;  r=l,...,R  the 
market  prices  of  all  resources  during  the 
period  t. 

(3)  Q*  =  ♦q(P1.  P2,  .....  PR). 

It  is  alos  assumed,  that  periodic  extractions  of  r  in  field 
j  are  subject  to  an  upper  bound  of  the  form. 

(4)  ijj  <  g£j(xrj>  Kp>    |f '  H  -  °   f°r  a11  r'  j  and  t- 
The  upper  bound  g  .  is  intended  to  reflect  the  impacts  of 
resource  and  capital  stocks  in  terms  of  impinging  on  periodic 
rates  of  extraction.   For  example,  given  a  capital  stock 
"K.",  (Pumps,  Pipe  size,  etc.),  the  smaller  the  recoverable 
stock  of  petroleum  in  petroleum  reservoir  (and  therefore,  the 
lower  may  be  the  pressure  for  natural  derive),  the  samller 
one  would  expect  the  maximal  periodic  rate  of  extraction. 

For  a  given  resource  stock,  smaller  capital  stocks  imply 
smaller  periodic  rates  of  extraction. 

Further,  it  is  assumed  that  r  extracted  from  j's  deposits 
are  used  for  increasing  stockpiles  of  r. 
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The  intertemporal  changes  in  stockpiles  of  the  resource 
r,  are  described  by  the  following: 

(5)  St+1  =  S*  .+W*  .+Mt  .-A*  ..,  S°.   . 

rj     rj   rj   rj   r j    rj  given,  for  all  r,  j  and  t. 

By  (5),  stockpiles  of  r,  in  field  j  at  the  end  of  period  t 

"S  .  "  equals  initial  stockpiles,  "S  ."  plus  extraction 
rj  rj  * 

"W  .",  plus  imports  "M  .",  minus  substractions  from  stock- 
piles "A  ."  for  current  consumptions. 

The  following  functional  forms  are  assumed  for  A  .  and 
&  rj 

rj 

(6)  A*  .  =  ^..(n*.'       T,  p*  .  .      T,  y*  .  .      T,  Q*), 
v     rj    r6jv  rj=l,...,J'  Krj=l, . . . ,J'  Mrj=l,...,J'  ^rJ' 

with  3^  >  0,  Jr   >  0,  ^  >  0,  ^  >  0 

(7)  Mrj  =  *mj(IIrj=l,...,J'  prj=l,...,J'  ^rj=l,...,J>  Qr> 

t  st 

where;  u    =  The  excess  of  P  ,  over  "C  ."  the  storage  cost  of 
rj  r'        rj  & 

resource  r  in  field  j,  during  period  t. 
Stockpiles  S  .  are  assumed  to  be  restricted  by  an  upper  bound 
N  .  which  reflects  the  capital  stocks  and  storage  facilities. 
With  K  .  defined  as  field's  J  capital  stock  for  storage  fa- 
cilities at  the  beginning  of  period  t,  it  is  defined. 

(8)  Sl .    <   N*  .(K* .),  ^i  >  0  for  all  r,  j  and  t . 

rj  -   rjv  rj"  3k  -  '  J 

In  the  resource  extraction  system  as  characterized  by 
(l)-(8),  two  observations  are  of  interest  at  this  juncture. 
First,  statement  (5)  allows  for  the  inclusion  of  port  facil- 
ities for  imports  to  be  treated  as  "resource  fields",  i.e., 
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*    t      t 
for  some  j=j  ,  X  .*,  W  .*  are  zero,  and  resource  management 

issues  in  terms  of  "field"  J   concerns  the  management  of 
stockpiles  of  r  from  imports.   Second  the  transition  equation 
for  stocks  given  in  (1)  is  in  a  simplified  form  for  exposi- 
tory purposes,  and  does  not  allow  for  the  dependence  of 
recoverable  stock  on  (earlier)  time-rate  of  production  and 
investment  which  may  be  of  particular  importance  for  re- 
sources such  as  petroleum.  [22] 

The  system  (l)-(8)  is  assumed  to  apply  to  fields  which 
are  known  (for  example,  "proven"  reserves  at  t=l)  at  the 
initial  time-period. 

To  allow  however,  for  exploration  expenditures,  e   , 

r  z 

which  may  result  in  new  fields,  z-l,...,Z,  with  deposits  of 

the  resource. 

Let  exploration  expenditures  e   result  in  discovery  of 

h   (e   )  units  of  the  resource  r.   (h,  is  essentially  a  meas- 
rz   rzy 

ure  of  the  "probably"  rewards  of  exploration  expenditures). 

By  assuming  a  concave  set  of  relationships  of  the  form: 

T      t 
(9)    I   hrz(erZ)  -  Xrz'  for  a11  r'  Z'  and 

9h        9h  ,    3h  0         9h  „ 

rz  >  o     rl  >    r2  >  rZ 

3e   —   '  9e  ..    3e  „    3e  „ 

rz        rl     r2  rZ 

By  (9),  initial  exploration  expenditures  are  for  h  1  (which 

yields  larger  "Probable  rewards"  per  unit  of  expenditure  than 

h   ,  z=2,...,Z),  until  a  bound  X  ..  is  effective,  after  which 
rz'     '    '  ' '  rl 

h  Q,...,h    are  applicable.   The  bounds  X    are  included  to 
it  &  r  Li  rz 

allow  use  of  a  continuous  model  as  well  as  to  allow  for 
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establishment  of  new  fields,  (with  reserves  X   )  as  a  result 

rz' 

of  exploration  expenditures.   It  is  also  assumed  e    to  have 
the  form. 

(10)  e*  =  4>  (n  ,  p  ,  n  ,  Q„), 
v  '        rz    e   rz   rz   rz   ^r 

3<J)  9(J)  9<J>  9<f> 

with  rf   >  0,    jf   <   0,    jf-   >  0,    and  -^   >  0. 

where;  II   =  expected  excess  of  market  price  of  r,  over  future 
extraction  costs  of  resource  r  in  field  z. 

p   =  expected  excess  of  market  price  of  r,  over  future 
import  costs  of  resource  r  in  field  z. 

n   =  expected  excess  of  market  price  of  r,  over 
capital  costs. 

Q   =  expected  demand  for  resource  r. 

A  set  of  equations  paralleling  (l)-(8)  is  thus  required 
for  "new  fields"  which  result  from  exploration  expenditures. 
The  relations  have  the  same  interpretations  as  (l)-(8)  but 
simply  apply  to  field  z=l,...,Z. 

(11)  Xt+1  =  Xt  -W*  +ht 

rz     rz   rz   rz 

<12>   Wrz  =  *wz<nrz<  ^rz'  Qr> 
(3>)   Q*   =  *q(Pl,  P2,..-.PR) 

t     t    t    t  for 

(13)  VT   <  gx  (X1  ,  Ku)  all 
v      rz  —  brzv  rz   z 

z      r 

(14)  St+1  =  St    +W*  +M*  -A1  ™d 
v   '    rz     rz   rz   rz   rz  t . 

(15)  Arz  =  *6z(IIrz=l,...Z'  prz=l,...Z'  ^rz=l,...Z'  Qr} 

(16)  MXrz   =    4>mz(n^z=1)  _    z,  PrZ=l,...Z'  *W=1,...Z'  Qr} 

(17)  S*  <    N1  (K*  ) 

rz  —   rz   rz' 
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Note  that  in  (11),  an  initial  stock  X    is  not  given  as  in 

rz 

(1).   If  expenditures  e   are  made  during  the  same  period  t, 

rz 

equation  (11)  becomes: 

(111)  Xt+1  =  h*  -W*   and  in  t+1 
v    '    rz     rz   rz 

(112)  xt+2  =  xt+1-wt+1+ht+1 

rz     rz   rz   rz 

and  a  "new  field"  z  is  brought  into  the  model. 

There  are  two  general  types  of  capital  stocks  which  are 

used  in  the  system.   Those  used  for  the  extraction  of  resources 

in  each  field  j    and  z,  K  ,  K  ,  and  those  stocks  used  for  the 

z        z 

management  and  storage  of  stockpiles  of  each  r  in  each  field 

— t    t  t      t 

j  and  z,  K  .,  K   .   In  each  field  j  and  z,  K.  and  K   are  lxN 
°  r j '   rz  J       '   j      z 

vectors,  where  K  .,  K   ,  n=l,2,...,N  are  intended  to  reflect, 

nj    nz 

first  a  quantitative  measure  of  various  specific  capital 
items  (pumps,  buildings,  machinery,  etc.)  and  second,  a 
qualitative  measure  of  capital  items  in  terms  of  technologies. 
Thus,  K_ .  may  be  a  petroleum  refinery  with  hydrocracking 

capabilities,  and  K   a  physical  plant  identical  to  K_ . 

'      nj    *   J  *  nj 

except  that  hydrocracking  capabilities  do  not  exist.   For  any 
n,  n,  marginal  investments,  substantively  different  from 
normal  capital-formation  notion  of  investment  (to  be  denoted 
as  V . ,  V  ,  later  in  this  model),  may  have  the  effect  of 

J     25 

"converting"  K-  .  capital  stocks  to  K  .  capital  stocks,  where 
&    nj    *  nj    * 

such  "conversions"  have  the  effect  of  "converting"  K_  .  capi- 

&    nj 

tal  stocks  to  K  .  capital  stocks,  where  such  "conversion" 
have  the  effect  of  changing  the  technological  structure  of 
the  extraction  processing  process. 
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Let  a_   .  be  the  investment  expenditures  which  converts 
n  nj 

K  .  capital  stocks  to  K_ .  stocks,  a   -.to  convert  K- .  to 
nj  nj  n  nj  nj 

K  . ,  a-    the  investment  expenditures  which  converts  K 
nj    n  nz  e  nz 

capital  stocks  to  K_   stocks,  and  a     -      to  convert  K-   to  K 

nz  n  nz  nz     nz 

capital  stocks.   Assuming  the  conversion  expenditures  to  have 
the  functional  form. 

(18)  o-       .    =   <b  -(Iu,  n*,  Q*) 

(19)  ax  - .   =  a  .(n*,  n*.  Qt) 

(20)  a-    =  *  -(nt,  nt,  Q*) 
v   '    n  nz    roz   z    z    r 

(21)  a*  -  =  4)  (n*,  n*,  Q*) 
'    n  nz   Tazv  z   z'   r' 

R  R 

where;  n*  =   E  (n*  ),   n*  =   E  (II*  ), 
J    r=l   rj     z   r=l   rz 

t    R   t      t    R   t 
nJ  =  r^(T1rj>>   ^z  =  ^/V^ 

If  f=cj).(a)  is  the  quantitative  measure  of  the  resulting  changes 
in  the  capital  stocks  after  conversion  expenditure  a,  then 

(22)  f-   .  =  d>„-r(0*   .) 

v      n  nj    yfjv  n  nj ' 

(23)  £*_.»*,  .(at  -  .) 
v   '    nnj    ^fjv  n  nj 

(24)  f-    =  d)--(a-    ) 

n  nz   yfz   n  nzy 

(25)  tX   -      =4.  (a1  -  ) 

n  nz    rfz   n  nz 

Define  V   as  additional  invest!   ;  to  capital  stocks 
which  are  used  for  the  extraction  of  resources  in  each  field, 
and  V   as  additional  investments  to  capital  stocks  which  are 
used  for  the  management,  import  and  storage  of  stockpiles, 
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and  having  the  general  formats: 

(26)  V*  .  =  <J>  .(II*,  p*,  y*   ni,  Q*) 
v  '        nj   Yvjv  j   J   J   J   r' 

(27)  vt  =  $  (n*,  p*.  y*,  n*.  Q*) 

v   '    nz   Yvzv  z'  Kz   z   z   r 

(28)  V*.  =  ^-.(n1.,  p*.,  y*   nt-,  Q*) 
v  J        rj   yvjv  rj'  ^rj   Mrj  '  'rj  ^v' 

(29)  v*  =  <f>-  (n*  ,  p*  ,  u*  ,  n*  ,  Q*) 

K      '        rz   Tvzv  rz'  Mrz'  Mrz'  'rz'  ^r' 

t    R   t 
Where ;   p  .  =  E  ( p  •. ) , 

J   r=l  rj 

t    R   t 
y  .  =   E  (y  .  )  , 

J   r=l  rj 

t    R   t 

t    R   t 

and   y  =  E  (y   ). 

r=l 
Also  define  in  general  format: 

(30)  D*  .  =  <f>  ..(W*  . ,  V*.,  K*  .) 

nj    Tdjv  r j    nj    nj 

(31)  D*   =  <j)  ,  (W*  ,  V*  ,  K*  ) 
'    nz   Tdz   rz   nz   nz 

(32)  D*  .  =  <f>,-(Vt  .,  K*  .) 
v  '        rj    ydjv  rj'   rj' 

(33)  D*   =  4>^-(Vt  ,  K*  ) 
v   '    rz   Ydzv  rz   rz' 

where  D  as  a  net  depreciation  function  for  extraction  capital 

gr\  gr\  grj  

stocks,  with  -tctt,   >  0,  -ktt  <  0,  ^s  >  0,  and  D  as  a  net  depreci- 

diV  —      dv  —      d  ii  — 

ation  function  of  capital  stocks  used  for  management  of 
stockpiles  with  — -  <_  0,    —  >  0-   Then  the  following  transi- 
tion  equations  are  assumed  to  describe  the  time  path  for 
capital  stocks: 
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N 
t+1     t    t       t      t 

(34)  IT.      =  K  .-D  .+  Z    (fx  _.-f_   .), 

nj     nj   nj  _=1v  n  nj   n  nj" 

nfn 

K  .  given,  for  all  n,  5=1, . . . ,N;  j=l,...,J;  t=l,...T. 

N 
t+1     t    t        t      t 

(35)  Kx  ■"■  =  Kx  -Dx  +  E  (fx  _  -fl    ), 
v      nz     nz   nz  -  -   n  nz   n  nz 

n=l 

for  all  n,  n=l,...,N;  j=l,...,J;  t=l,...,T. 

(36)  K^1   =   K^j-D*..,  for  all  j=l,...,J;  t=l,...,T, 

Kr .  given 

(37)  Kt+1  =  K*  -D*  ,  for  all  z=l,...,Z;  t=l,...,T. 

v     /      rz         rz    rz>  >  >        >  >  » 

A  word  of  explanation  may  be  required  at  this  point  in 
order  to  clarify  the  meaning  of  (34)-(37)  above.   The  first 
two  of  the  transition  equations  specify  that  the  capital 
stock  of  type  n  available  for  use  in  extraction  field  j 
(or  z)  at  any  point  in  time  depends  on  the  size  of  the  stock 
in  preceding  period  less  net  depreciation  over  the  preceding 
period  (which  may  actually  be  negative  if  investment  in  new 
stocks  is  large  and/or  if  extractions  are  small),  plus  the 
net  conversion  of  other  types  of  capital  stocks  to  type  n. 
Thus,  the  latter  summation  term  takes  account  of  the  amounts 
of  other  classes  of  extractive  capital  (n^n,  h=l,...,N)  con- 
verted to  type  n,  less  that  of  type  n,  converted  to  other 
types. 

Equations  (36)  and  (37)  are  required  for  the  available 
capital  used  to  manage  and  store  the,  up  to  R,  different 
energy  resources  in  field  j  (or  z)  at  time  t+1.   Since  it  is 
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essentially  assumed  homogeneity  of  this  class  of  capital  for 
any  given  type  of  energy  resource,  there  are  no  conversions 
to  quantify. 

Next  define  U   as  the  quantity  of  energy  resource  r 
available  for  use  in  the  production  of  "final  energy  commod- 
ities" at  time  t. 

Thus  U  may  be  thought  of  as  the  volume  of  petroleum 

(crude)  made  available  to  refineries,  coal  at  trackside,  etc. 

The  term  U   then  includes  the  sum  of  draw  downs  from 
r 

stockpiles  of  r  from  all  resource  fields,  j,  and  z,  for  cur- 
rent consumption,  then; 

t    J   t   Z   t 

(38)  U   =   I  Au.+  £  Au  ,  for  all  r=l,...,R;  t=l,...,T. 

r    j=l  rj  z=l  rz 
Next  it  is  assumed,  there  exist  L  different  "final  energy 
commodities"  and  define  A.,  as  the  amount  of  the  1    type 
produced  at  time  t,  (1=1,..., L).   The  amount  of  each  produced 
depends  upon  the  amount  of  the  R  possible  energy  resources 

used  as  inputs  "U   ",  as  well  as  up  to  Q  different  kinds  of 

~t 
capital  stocks  "K  ,"  used  in  the  production  process.   Thus 

the  production  functions  could  be  written  as: 

(39)  A*  =  ♦1(uJ1....fU*1;  K^,...,^),  for  all  1=1,..., L; 

t=l      T 
9A.T       BA^ 
with  — £  >  0;  — q:  >  0  for  all  r,  q. 

Note  that  obviously  all  energy  resources  might  not  be  used  as 

8A 
inputs  in  each  production  function  process;  then  -r= —  may  be 

3Url 

zero  for  some  r  in  each  1. 


43 


Now  let  each  unit  of  energy  commodity  1  contain  a  BTU' s . 
If  E   equals  the  minimum  amount  of  BTU  required  in  period  t, 
then  the  problem  below  must  be  constrained  by: 

L    t     t 

(40)  Z      ,A   >  E  ,  for  all  t=l,...,T. 

1=1  X    L 
Let  H,  j>  0  equal  the  minimum  required  quantity  of  the  1 
specific  energy  commodity  in  period  t.   This  implies. 

(41)  A*  >  H*,  for  all  1=1, . . . ,L;  t=l, . . . ,T. 


C.   THE  CRITERION  FUNCTION 

The  energy  independence  policy  presented  by  President 
Nixon  has  been  usee1,  as  an  application  of  the  model  introduced 
in  previous  section  and  for  the  formulation  of  its  criterion 
function. 

The  basic  concept  of  this  policy  arises  from  the  threats 
to  the  security  of  the  nation  because  of  foreign  supply 
interruptions,  and  the  assumption  that  the  threats  most  likely 
to  be  encountered,  are  best  met  through  domestic  self-suffi- 
ciency.  Therefore,  it  calls  for  more  storage  of  resources, 
their  higher  production,  and  expansion  of  explorations  for 
new  deposits  and  sources  of  energy,  aimed  toward  the  realiza- 
tion of  this  concept  by  period  T.  (1980)   Within  the  decision 
situation  described  by  the  conditions  and  restrictions  (1)- 
(41)  of  preceding  section,  concern  here  is  with  the  indirect 
selection  of  values  for  extraction  rates  of  natural  resources, 
the  expenditure  rate  for  exploration,  rate  of  investment  of 
capital  stocks,  rate  of  production  of  energy  commodities  to 
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respond  to  energy  requirements  and  production  such  that  the 
following  expression,  the  present  value  of  total  capital  out- 
lays and  consumer  costs  over  "T"  periods  of  time  be  minimized, 

T 

(42)  Min   E  C* Bt+i|;T(XT,KT)3t 

t=l      X 

Subject  to  (1)-(41) 
where;  C   =  Sum  of  all  capital  outlays  and  consumer  costs. 
ty~   =  A  concave  terminal  value  function  for  all  the 

natural  resources  and  capital  stocks.  [23] 
3   =  l/(l+r),  where  r  is  the  appropriate  positive 

discount  rate.   Reflects  the  discounting  factor. 
The  objective  function  in  a  detailed  form  could  be  written  as 

T 

(43)  Min   E  {E  et  Bt+S.  V*  .6t+Z  V*  8t+E.Vt.6t+Z  V*  B* 
v   '       ...  rz  rz   nj  nj   nz  nz   rj  rj   rz  rz 

+E  oX   _.Bt+E  a1  -  Bt+E.  a*  .^+1   o-        Bt 
-.  n  nj   _   n  nz   nj  n  nj   nz  n  nz 
ii  j         nz 

+E.  W*  .•Cwt6t+E  W*  *Cwt  Bt+l.  Sl  .•CStBt 
rj  rj   rj   rz  rz   rz    rj  rj   rj 

+  £  St  •CStBt  +  E.  Mt .•CmtBt  +  E  M*  »Cmt&t} 
rz  rz   rz   rj  rj   rj   rz  rz   rz 

T   T   T 
+(J;T(X\K1)B1. 

Subject  to  equations  (1)-(41) 
where ; 


T      T     T 
X   =  E.  X  .  +  E  X 

rj  rj  rz  rz 

K   =  EK  .  +  EK  +E.  K1  .  +  E  K   +E  K,  . 
j  j  z  z  rj  rj  rz  rz  iq  lq 


R 

Z 

E      • 
rz 

=       E 
r=l 

E   ; 
z=l 

R 

J 

N 

J 

E.     ■ 

=      E 

E   ; 

E     ■ 

=      E 

E 

rj 

r=l 

j=i 

nj 

n=l 

3=1 
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D.   EQUILIBRIUM  CONDITIONS 

The  optimization  problem  expressed  by  (43)  and  equations 
(1)  through  (41)  could  be  solved  by  Lagrangian  method.   As- 
suming that  the  functions  in  (1)  through  (43)  are  continuous 
and  possess  first  partial  derivatives  with  respect  to  each 
argument,  the  following  Lagrangian  expression  may  be  formed. 

T 

(44)   L  =  E  {E  e*  Bt+Z.  Vt.3t+E  V*  Bt  +  Z.  V*  .  B^E  V*  6t 
v  '  .  -  rz  rz   nj  nj   nz  nz   rj  rj   rz  rz 

tttttttt 
+  E   a1   _.B   +Z   ax   _    B   +E.  a-      .3  +1   a-        B 
_  .     n  nj       _.     n   nz       nj     n   nj       nz    n   nzK 

+  Z.  W*  .•Cwt6t  +  E  W*  •CWtBt  +  E.  St  .•CSt6t  +  Z  S1  •CStgt 
rj  rj   rj   rz  rz   rz   rj  rj   rj   rz  rz   rz 

+  Z.  M*  .•Cmt6t  +  E  M*  •CmtBt-1At  .  [X*  .-X°  .+  £  f*.] 
rj  rj   rjM  rz  rz   rzp    r j  L  rj   rj    1  r j J 

-E.  2At  .  [W*  .-(J)1  .(n*  .,  p*.,  Q1)] 
rj   r j l  rj  vwj   r j   Krj   ^ryj 

-Z3Xt[Qt-c()t(P,  ,  P,  , P~)] 

r  r  r  q   i    2'     Q 

-Z.   "X*  .  [W*  .-gt  .(X*  .,  K1)] 
rj   r j   r j  6r j   r j    j  J 

ct    t    o    •*   +      t    *   t 
-5A.[S  .-S  .-  E  W  .-  E  ML.+  E  A  .] 
rj1  rj   rj  t=1  rj  x=1  rj  x=1  r j 

-E.  6Xt  .  [A*  .-d>*  .(n1  .  i     ..p*.  i     T,y*   i    t.Q*)] 
rj   r j l  rj  T6jv  rj=l, . . . J' Krj»l, . . . J' Hrj=l , . . J'  "r/ J 

-e.  7xt  .[M1  .-4)*  .(n* .  ,   t.p*.  -   t.u*.  1   T,Q*)] 

rj   r j L  rj  ymjv  rj=l, . . J' Hrj=l, . . J' Hrj=l, . . J'^r  J 

-  E.  8Xt  .  [S*  .-N*  .(I*  ■)]-9Xt    [ht  (e*  )-X   ] 
rj   rj   rj   r j v  r j     rz1  rzv  rz7   rzJ 

-E  ^X*  [e*  -0  (nt  ,q1    ,r\t    ,Qt)]-6t 
rz    rz   rz  Ye   rz  rz  'rz  r 

-  llX*  [X1  +  E  W1  -  E  h1  ] 

rz   rz    i  rz    -  rz 

T  =  l      T=l 

1   2   i   t        rTirt  ±  /rrt  _t  /~vt 


•E  ,ZX"  [W  -4)   (ir  ,  p"  ,Q  )] 
rz    rz   rz  wz   rz   rz   r 

■  E  l  3Xt  [W*  -gt    (Xt  ,  K*)] 
rz    rz   rz   rz   rz    z 
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-lkXt    [S*  -  E  W*  -  E  M*  +  E  A*  ] 

rz   rz    -  rz    -.  rz    -  rzJ 

T  =  l      T=l      T=l 

-E  15Xt    [A1  -<|>*  (n1   ,    ..p*   ,    -.y*   .    -.Q*)] 

rz    rz  rz  v6z   rz=l,,.Z  rz=l,..Z'  rz=l,..Z'  rn 

-E  1BX*  [M*  -(J**  (n*   -    „,pt   .,    -.y*   -    -.Q*)] 
rz    rz1  rz  Ymz   rz=l, . . Z' Krz=l, . . Z' ^rz=l , . . Z'^r' J 


-E  I7X*  [S*  -N*  (Kt  )]-Z.  ^A*  .Bt[a^   .-<j>  ,(n*  .n1  .Q* )  ] 


-I  19X*  [a*  _.-*  .(n*  n^,Ql)]6t 


nj 


njl  n  nj  Taj^  j   j  ^v 


-E  20Xt    [a-   -<t>  ^(nt,nt,Qt)]6t 

nz    nzL  n  nz  TaZv  z'  'z'  r/J 

-E  21X*  [a*  _  -(f)   (nt,nt,Qt)]6t 

nzL  n  nz  Tazv  z' 'z'  r/J 
nz 

-E22X*[f-   .-()>.- (a-   .)]-E  23X*[f*  -.-<{>.,  (a*  _.)] 


nj    n j   n  nj  Yfjv  n  nj 


nj 


nj   n  nj  Tfjv  n  nj 


JL  2*X*  [f*   -<|>*=(aj;  „„)]-£  25xLff!  s„-*^„(a*  -  )] 


nz    nz   n  nz  Tfz  n  nz 


nz 


lnzl  n  nz  Tfz   n  nz' 


-E.  26Xt  .Bt[Vt  .-<J)  .(n^p^.y^.Ti^.Q*)] 
nj    n;T    nj  Yvjv  j'kj,hj'  'j'^r/J 

-j  27X*  Bttvt  -<j)  (nt,Pt,vit,nt,Qt)] 

nz    nz   L  nz  rvzv  z  Kz   z   z   z/J 

-E.  28Xt.gt[Vt  .-())_  .(n*  .,pt  ..y*  ..n1  ..Q*)] 
rj    rjH    rj  rvjv  r  j  "r  j  '  pr  j  '  'rj'^r'J 

-E  29xt  3t[vt  -<j)-  (n*  .p*  ,yt  .n*  .Q*)] 

rz    rz   l  rz  Yvzv  rz  ^rz'  rz'  'rz'  ryj 

-E.  30Xt  .  [D1  .-4)  ,  .(W*  .,  V*.,  K*.)] 
nj    nj   nj  ydjv  r j '   nj '   nj  J 

-E  31X*  [Dt  -<J>,  (W*  ,  V*  ,  K*  )] 
nz    nz   nz  Tdz   rz'   nz   nz  J 

-E.  ,,At.[Dt  .-A^-fV*  ..K*  .)] 

rj    r j   rj  ydjv  r j   r j 

-E  ''X*  [D*  -(K-CV1-  .K*  )] 
rz    rz   rz   dz   rz   rz 

-•*X*  4[K*  .-K°  .-  E  D1.-  E    E  (ft  -A-t-       .)) 

nj1  nj   nj  t=1  nj  t=1  -=1^    n  nj   n  nj 

t       t    N 

-35Al  [K   -  E  D   -  E    E  (fL  _  -f_    )] 
nz   nz   1  nz   ..  _  - v  n  nz  n  nz 
t=1     t=1  n=l 
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-36Xt.[Kt.-K°  .+  Z   Dt.]-37Xt    [K1    +  Z  D*  ] } 
rjl  rj   rj  t=1  r j     . rzL  rz  x=1  rzJJ 

+j¥{-38AltAr*l(Url'---URl'  ^r---^l>] 

-39A*[  Z  a1A*-Et]-40A^[A^-H*]}+  ^T(XT,KT)3t. 

The  necessary  conditions  for  minimizing  equation  (44) 

include  the  following  equations,  which  are  then  respectively 

interpreted  for  their  economic  implications  on  the  objective 

function. 

3(Jt 

(45)  -Q-  =    6t-26Xt.6t-30Xt.[-4i] 
3Vt  nj       njl3Vt  J 

nj  nj 

In  equation  (45)  the  term  26A  .  6   represents  the  discounted 
value  of  the  incremental  change  in  objective  function  due  to 
a  unit  of  change  in  the  investment  of  capital  stocks  used 
for  extraction  of  resource  r  in  field  j,  defined  by  Constrain- 
ts   ^H  ' 

ing  equation  (26).   The  terms  30A  .[ — r**-]  represents  the 

nj  ayt 

nj 
implicit  value  of  a  marginal  change  in  net  depreciation  of 

these  capital  stocks  due  to  a  marginal  change  in  relative  in- 
vestments.  The  necessary  condition  for  optimality  requires 

that  -~-  =  0,  or 
nj 

(45-1)  30\t.[ — P-]    =    pt-a6Xt-Bt.   ' 

nj1   t  njp 

nj 

A  similar  interpretation  can  be  made  for  equations  (46)- 
(52)  below. 
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(46) 


9L 


3V 


=  6t-27At  B*-31**  [—p]   =   0 


nz 


nz 


nz  3yt 
nz 


86 

(46-1)   31Xt  [—Is]  =  B*-27**  B* 

v     '      nz  3„t  nz 

nz 


(47) 


3L 


=  Bt-*«XtJ3t-8»Xt.[-43L]  =  0 


3V. 


rj 


rj 


rJ  3V*. 


84  .  -t 


(47-1)   32X^t^]  =  B*-2^.  B* 

rj 

»\  j.  ** 

(48)   -4t-  =  B*-2^  Bt-33Xt  [-^]  =  0 

3V  rZ       rz  3V 

rz  rz 

(48-1)   33Xt  [— |i]  =  6t-29Xt  B* 

rz  3V1  rz 

rz 


(49) 


-^ =  Bt-19Xt.Bt  +  23A^.[  P       ] 

3at  -.         n^      "J^  - 


n  nj 


n  nj 


36*  . 


(49-1)   23X*.[ — p — ]  -  19Xt.Bt-Bt 

n  nj 


(50) 


3L 


30*  - 
n  nz 


=  0 


t   ^■f 
(50-1)   25X-  [ ^ — ] 

da  - 
n  nz 


B*!2^-  -1] 
nz   J 


(51) 


3L 


3a!   . 
n  nj 


=  0 


t 

36  --: 


(51-1)   22Xr  .[ JJ— 

nJ  3a^   . 
n  nj 


6t[18Xt .-1] 
nj   J 
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(52)  — ^—  -  0 

n  nz        . 

(52-1)   "X*  [ — P — ]  =  Bt[20Xt  -1] 
v     '     nzl~  t    J      L    nz   J 
3a_ 
n  nz 

gxt 

(53)  -^-  =  Bt-Cw*-1Xt.-*Xt4-*Xt.+8Xt.-80Xt.[— P-]    =  0 

v   '   3wt        rj    rj    rj    rj    rj    nJ  9Wt  J 

rj  rj 

t  wt 
In  equation  (53)  the  term  B  «C  .  represents  the  discounted 

t   ^di 
value  of  the  extraction  cost,  and  30X  .[ — -zr-]    represents  the 

nj  9Wt 

rj 

implicit  value  of  the  marginal  change  in  net  depreciation  of 

the  capital  stocks  due  to  a  marginal  change  in  extraction 

rate  of  resource  r.  where  *X   represents  the  value  of  the 

rj    * 

incremental  change  in  objective  function  due  to  a  unit  of 

change  in  the  extraction  rate   restricted  by  equation  (1) 

and  2X    for  equation  (2)  and  1*X    for  equation  (4)  and  5X 

rj      *  rj      M  rj 

for  equation  (5).   The  necessary  condition  for  optimality 
requires  that       . 


(53-1)   3  0Xt  f   dij  =  Bt.cwt-(ixt.  +  2Xt.  +  nt.-5Xt.) 
njL3wt         rj    rj    rj    rj    r j 

rj 
A  similar  interpretation  could  be  made  for  the  equation 
(54)  below. 

(54)  -Q-  =    Bt'Cwt-11Xt   -12X*  -iaX*  +  lkXt   -3lXl    [— P-]    =   0 
v   '   awt        rz    rz    rz    rz    rz    rzL~wt  J 

rz  rz 

(54-1)  3lXl    [—4^]    =    Bt-CWt-(11Xt  +12At  +  13At  -^X*  ) 
v     '  rz  „T„t         rz     rz    rz    rz    rz 

d  W 

rz 

9ht 

(55)  -Q-  =    Bt-9Xt    [-45]-10At    Bt  +  11Xt      =   0 
~    t  rzL.tJ  rz  rz 

ae  3e 

rz  rz 
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In  equation  (55)  the  term  3X      [ — — ]  represents  the  implicit 

rz  r>.  l 

rz 
value  of  a  marginal  change  in  units  discovered  of  the  re- 
source r  in  field  z  due  to  a  marginal  change  in  the  explora- 
tion expenditure  restricted  by  equation  (9).   The  term 
10  A   3   represents  the  discounted  value  of  the  incremental 
change  in  objective  function  due  to  a  unit  change  in  explora- 
tion expenditure  described  by  equation  (10).   Similarly  the 
term  11X        represents  the  value  of  an  incremental  change  in 
objective  function  due  to  a  unit  change  in  the  quantity  dis- 
covered of  resource  r  described  by  equation  (11).   The 
necessary  condition  for  optimality  requires  that 

.   ah        ,       .       . 

(55-1)   9Xt  [—£2]  =  Bt(l-10At  )  +  11Xt 

v     '    rz*   t  J  v     rzy     rz 
3e 
rz 

The  most  important  consideration  in  this  study,  is  the  impact 
of  price  changes,  on  the  objective  function  and  other  vari- 
ables in  the  model.   This  is  described  by  (56),  which  is  the 
partial  derivative  of  equation  (44)  with  respect  to  P  ,  the 
price  of  resource  r  in  period  t, 

„  ,     a^.  an*.      a$*     ap*.        .   a^ 

(56)    i\  =  .x*[-!=i.-|i  +  -!EL.^|i]  +  .xJI-!a] 

apL         rj  an  .  ap         ap  .  ap  ap 

r  rj         r  Krj         r  r 

.     3<j>   .  an  .      a<j>..  ap   .      a<t>..  au   . 
+&x     r 2J. Ll  +  — °I LI  +  — °1 Hi 

ri  l       t  t  t  t  t  t    J 

3  an*..  dPl      3Pv  a?*      dV;.  dPl 

t       3*mi    3IIri         3*mi    3pri         3*mi    3vJri 
+  7Xt    i — HL1 Ll  +  — Uli Ll  +  — EI Lli 

rilt  t  t  t  t  tJ 

J  an   .  ap  ap    .  ap  au   .  ap 

rj        r  Hrj        r  Mrj        r 
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10X*_[ 


d^ 


rz 


12X*_[ 


3<j> 


wz 


r'l3nt 

rz 


15xL[ 


3$ 


6z 


rz 


+  16xL_[ 


84) 


mz 


'-'an* 

rz 


an 


rz 


3P 

an 


r 

t 

rz 


3P 

an 


t 

r 

t 
rz 


3P 

an 


t 

r 

t 
rz 


3P 


3<J)         3p 


3p 

3<t> 


rz 

t 
wz 


3p 


3<J> 


rz 

t 

6z 


3p 

3<j> 


t 
rz 

t 

mz 


rz 


3P 
9p 


r 

t 
rz 


3P 
3p 


r 

t 
rz 


3P 
3p 


r 

t 
rz 


3p~      3p 


3*       3n 


3n 


rz 


ay 

3<J> 


rz 

t 
mz 


rz 


ay 


3P" 


rZ]3* 


34,6z    3y 


rz 


3P 

ay 


r 

t 
rz 


rz 


3P 


t      t    ScK-    311    , 

J      an;.  3p; 

rj        r 


9Vj    9Pr 


t     t  3(f).,  sn    . 

+  "xSjB    t— P — P 

J       an   .  3P 

rj        r 


9*6j_9nrj 

anjj  ap* 


.2  0 


.       .     34.-    3n 

t    _t r    r6z 


*nz8    1 


2,4eti 


3n 


3<J) 


rz 


6z 


sn 


■26X^B"[ 


rz 


.     3*     . 
t  r    vv^ 


nj 


an. 


rj 


rz 


3P 

3n 


r 

t 
"rz 


3P 

sn 


t 

r 

t 


aP 


Y6z        rz. 


3n 


rz 


3<J) 


6z 


an 


rz 


3<f> 


vj 


ap 

3n 


r 

t 
rz 


3P 

ap 


t 

r 

t 
rj 


3p     .    3P 
wrj         ] 


^llm^£l 


ay 


rj 


ap 


34,       an 

-¥■— v+ 

3n     .    3P 
'rj         r 


3d)         3u 
vvz      Mrz 

■--■■-    ■  •  

3yL  *pl 

rz        r 


r  ,  .  2  7 , t    „t .  "Tvz    ""rz 


a*, 


an. 


3d)-  •     3p     ■ 

3P1.  aPt 

prj    °    r 


Anz6    t 


9*vz    9Prz 


an1     3P1 

rz        r 


3pL    3P^ 
Krz        r 


!!vz.!^z       t  6tr!^i.!4j 

t    "      t    J+      ArjB    L„t      .„t 


an"     a? 

'rz        r 


t  t 

3d>-  •     an     .  3d)-  .     3n     • 


3nt.     3P* 
Hrj         r 


3n*-    3P1 
°"rj    °   r 


9nrj    9Pr 
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.        .     3<t>-       3H  3d>-       3p*  3d)-       3y* 

^2  9-vt    .t.    Tvz        rz    ,       Yvz      Hrz    ,       Yvz      Hrz 

rz    3n*   3P*     3p*   3P*     3u*   3P* 
rz   r      rz   r      rz   r 

3(j>-  an* 
+  _!^._^  =  o 

3<4  9Pr 
And  after  the  arrangement  of  the  terms,  it  could  be  written 

3n  .       3d)  .      .   3d),.  .      .   3d) 

(56-1)  -^1[2A*.^  *  U*  •-%  *  7A^-^ 
apC     rj    x       rj    x       rj  3nx 

r         rj  rj  rj 

,   ,  3d)  —       j.   j.  3d)  •     3n„      .   3d) 

+„xt  ,t.Jai+  ,,xt  (t._!2i   ♦_£s[».,£ z~!f- 

nJ       snr .  nJ       an.        3Pt  rz  3nr 

rj  rj  r  rz 

+  12    t    ^wz   +    l5    t    .^6z   +    16    t      ^mz   +    20    t      ^oz 
"•"J.  "<  "<  "^"rz 

+*,xt    t .*h  +  -^et.Hvz  +  „xt    t.^ 

3"rz  3"rz  »rz 

3P1    '      rJ    3P1.  rJ    »,,*.  «B      »„*  . 

r  Mrj  Mrj  prj 

t      t    3(f)-,  3nri  t      34>p  t      t    8<})rM 

+  28X*Bt.-^]    +-fL[1°X*    — 1-+    ,8X*B  *— f- 

3        »v5j  <  '    3t4  nJ        3n*j 

t      t    8*rri  t      t    3*vi  t      t    3<t>vi 

3nrj  3nrj  3nrj 
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3p                 ,       3d>  4.3<))  3d> . 

+  rz[10>t e_  +    12xt    .  -    wz   +    !5Xt    .    T6z 

~Dt  rz   ~    t  rz    ~    t  rz   ~    t 

dP  dp  dp  dp 

r  rz  rz  rz 

.       3<{>  .  3d>  ,  3d)  — 

+16xt   mZ  +  2?xt  gt._^VZ  +  29xt  gt.J^VZj 

rz  3pt       nz   apt       rz   apt 
rz  rz  rz 

By*      .   3t()t  8(f)  t  9d) 

+  I^r^xt  .   6z  +  iext  .  vmz  +  2?xt  pt.  vvz 

0T)t     rz  ~  t       rz  ~  t       rz    ~  t 
3P  3u  3u  3u 

r  rz  rz  rz 

.   ,  3d)-      3n  3d)  .  3d)  — 

+  29xt      gt.JVZ]      +    _^[!0Xt      ._!|_    +     20Xt      ^ OZ 

rz    «>.  l      ^-r»*t      rz  r.  x       rz    ^^t 
3y       3P  3n  op 

rz      r  rz  rz 

.   ,  3d)         .   ,  3d)  34- 

+  21Xt   3t._^OZ  +  27xt    t VZ  +  29xt    t VZj   = 

nz     t       nz     t       rz    a  t 
rz  rz  rz 

.  3<})t 

r  gpt 
r 

In  equation  (56-1)  the  terms  outside  of  the  brackets  repre- 
sent the  marginal  change  in  the  excess  of  markets  price  of 
resource  r  over  different  costs  in  different  fields  due  to 
a  marginal  change  in  the  market  price  of  resource  r.   And 
the  terms  inside  the  brackets  (some  of  them  discounted) 
represent  the  implicit  value  of  the  marginal  changes,  in 
different  variables,  influenced  by  the  price  of  resource  r, 
due  to  a  marginal  change  in  the  excess  price  over  the  costs 

In  optimality  the  total  sum  of  these  terms  in  the  left 
hand  side  of  equation  (56-1)  must  be  equal  to  the  term  on 

.  3d)t 
the  right  hand  side  -  A  [ — f]    which  represents  the  implicit 

r  3Pr 
r 

value  of  the  marginal  change  in  units  demanded  of  resource 
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r,  due  to  a  marginal  change  in  its  market  price,  the  demand 
function  is  described  by  equation  (3). 

The  above  mentioned  conditions  simultaneously  solved 
with  the  following  equations  would  determine  the  values  of 
multipliers  and  the  optimum  rates  of  exploration,  investment, 
extraction,  importation  and  conversion  expenditure. 

(57)  _ik_  =  ^.cmt.    +  5Xt.  -  7Xt.  =  0 

9M*        rj     rJ     rj 
rj 

(57-1)   5Aj\  -  7A^  =  -&*•<£* 

(58)  _lk_  =  ^.Cmt   +  ^X*  -  16X*  =  0 
3Mt        rz      rz      rz 

rz 

(58-1)   16X*   -  "X*   =  B^C™* 
v     '  rz      rz       rz 

(59)  -~-  =  ^'CSt   -    5Xt:    -  8Xt.  =  0 
3gt        rj     rj     rj 

(59-1)  5xt .  +  8xt .  =  s^c31: 

'    rj     rj       rj 

(60)  _ik_  =  Bt-Cst  -  "X*  -  I7X*  =  0 
v   '   3St        rz      rz      rz 

rz 

(60-1)   ^X*   +  "X*   =  ^-CSt 
v     '     rz      rz       rz 

(61)  |y  =  0;  for  all  Xs. 
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V.   CONCLUSIONS 

An  attempt  has  been  made  to  bring  together  and  formalize 
some  of  the  factors  that  seem  to  exist  in  "Supply  and  Demand 
for  Energy".   After  a  brief  review  of  the  world  energy  con- 
sumption it  is  shown  that  there  is  a  broad  parallel  between 
economic  growth  and  energy  consumption.   It  is  important  to 
explore  the  relationship  as  precisely  as  possible  and  thus 
more  intensive  research  is  needed  to  determine: 

a.  The  impact  on  energy  consumption  of  technological 
efficiencies  in  consuming  countries. 

b.  The  "average"  nation-wide  levels  of  energy  consump- 
tion required  for  given  level  of  GNP. 

c.  The  effect  of  above  mentioned  topics  studied  in  terms 
of  specific  energy  sources.  In  addition,  the  substi- 
tutability  of  one  energy  source  for  another  in  satis- 
fying requirement  should  be  considered. 

d.  A  more  specific  effect  of  relative  prices  of  resources 
on  energy  consumption. 

e.  Impacts  of  environmental  controls  and  pollution 
abatement  upon  the  volume  of  energy  demanded  and  upon 
shifts  to  substitute  energy  sources. 

f.  Impacts  of  economic  and  financial  policies  such  as 
balance  of  payments  or  balance  of  trade  considera- 
tions upon  the  volume  of  energy  consumption. 

By  introducing  the  above  mentioned  considerations  in  the 
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empirical  analysis  of  past  data,  a  greater  insight  might  be 
gained  into  future  projection  of  energy  requirements. 

In  the  analytical  section  of  this  study  the  Cummings  and 
Whipple  model  was  extended  by  posting  a  relationship  between 
extraction,  storage,  import,  exploration  and  market  prices 
of  resources  and  relative  cost  factors. 

The  model  does  not  include  social  values  and  environmen- 
tal considerations. 

The  model  would  allow  for  price  regulatory  policy  deci- 
sions, and  considers  both  supplies  and  demands  of  resources, 
and  their  interactions  is  achieved  through  the  dynamic  mar- 
ket mechanism.   It  is  designed  for  repeated  use  and  consid- 
ers the  technological  impact  and  conversion  of  capital  stocks. 

The  model  is  disaggregated  by  regions,  covers  all  energy 
forms  and  takes  into  account  the  interfuel  substitution 
brought  about  by  price  changes. 

It  is  believed  that  expansion  of  the  model  to  cover  the 
following  additional  areas,  would  make  it  more  useful  in  con- 
junctive management. 

1)  Inclusion  of  environmental  considerations  and  pollu- 
tion abatement. 

2)  Inclusion  of  specific  functional  forms  for  extraction, 
import,  investment,  etc.   These  could  be  obtained  by  statis- 
tic estimation. 

3)  The  inclusion  of  uncertainty  in  exploration  expendi- 
tures. 

4)  Inclusion  of  transportation  problems. 
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5)   Refinement  of  the  model  so  that  a  computer  simulation 
might  be  employed. 
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APPENDIX  A 

1.   Compound  interest  and  rate  of  growth 

If  $y0  is  invested  at  compound  interest  at  r%  per  year, 
compounded  yearly  then  the  amount  after  t  years  is  $y  where: 

(1)  y  =  Fo(l  +  r)* 

if  the  interest  is  added  n  times  a  year,  then 

r  nt 

(2)  y  =  y0(l  +  f) 

Discontinuity  is  an  essential  feature  of  this  compound  inter- 
est problem.   The  more  frequently  is  interest  added,  the 
larger  is  the  amount  of  a  given  sum  at  the  end  of  any  period, 

Letting  n,  the  number  of  the  frequency  of  compounding 
interest  to  take  larger  and  larger  values,  it  is  clear  that 

1  n 

(1  +  — )   tends  to  a  definite  limit  and  it  can  be  shown  that 

1  n 
e  =  lira  (1  +  -)  ;  e  =  2.71828 
_       n 

n->-o° 

Returning  to  the  case  of  compound  interest 

nt  n  rt  1    rt 

Y  =  yod  +  f )    =  yoUl  +  ^)r}    -  y„{(l  +  £)  } 

1  m 
for  n^  oo  ;  (1  +  — )   =  e   then 
v    m 

(3)  y  =  y0ert 

We  have  now  derived  a  concept  for  continuous  interest. 
Note:   This  model  represents  a  constant  percent  rate  of 
growth  (interest  rate).   It  is  clear  that  the  percent  rate 
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of  growth  at  time  t  is  r(t)  =  — — (100)  where  -?*r  represents 
the  rate  of  change. 

2.   Present  Values 

a.   If  the  interest  is  added  once  yearly  at  r%  then  y 
is  the  present  value  of  y  available  after  t  years;  from  Eq. 
(1): 

(4)     y0  =  — Z-- 


(l+r)t 


b.   If  interest  is  added  n  times  a  year  at  r%  from  Eq 
(2): 

(5)     y0  =  —  " 


r  nt 

(l+£) 
v   n 


c.   Finally,  if  interest  is  added  continuously  at  r% 
from  Eq.  (3)  it  implies 

(6)     y0  =  -%p£     or  y0  =  ye~r 
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APPENDIX  B 

Calculus  of  variations  -  Euler's  Equation. 
If  the  extreme  values  of  the  integral 

u  =  /  f(t,q(t),|a)dt 

are  required  for  all  possible  variation  in  the  function 
Q=q(t)  such  that  q(t0)=qo  and  q(t1)=q1  where  (t0,q0)  and 
(ti.qi)  are  fixed  points.   The  function 

f(t)  e  f(t,q,||)  e  f{t,q(t),q(t)} 

which  gives  variable  u  on  integration,  depends  on  the  variable 
t,  on  the  variable  function  q=q(t)  and  on  the  first  deriva- 
tive q(t).   This  is  a  function  of  t  given  in  the  function  of 
functions  form.   In  solving  the  problems,  the  function  q(t) 
is  taking  in  the  form  q=q(t ; a, 6, y> • • • )  where  q(t ; a, B, y, • • • ) 
is  assumed  to  be  a  fixed  function  (with  a  continuous  deriva- 
tive) and  where  a,  $,y , . . .    are  parameters.   Given  differential 
increments  6a,  66,  6y>  ...  to  the  parameters,  we  derive  the 
corresponding  variations  5q  and  6q  in  the  function  q  and  its 

derivative  q  =  -rr    '• 

dt 

<*>         si  =  £  6«  +  !t  «  +  |a6lf  ♦  ....  and 

■  &<£>«« ♦  SH*>«  ♦  &<f>  + . . . 


1 


Reference   7,    Page    521. 
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The  symbol  "6"  is  used  in  order  to  distinguish  parametric 
differentials  from  "d"  deferring  to  variation  in  variable  t. 
The  function  f(t,q,-r^)  and  the  integral  u  can  now  be  consid- 
ered as  dependent  on  the  parameters  a,  &,y, . . .    and  the  varia- 
tion in  their  values  are  obtained  as: 

from  Eq.  (2)  results 

(3)      8f=|§«q  ♦■§&(*> 

and 

ti       .  t, 

6u  =    /  f(t,q,q)dt  =  ./  6f  dt 

to  to 


From  Eq.  (3) 


t]  ^        t 


(4)  6U  =  J  (||6q)dt  +  /{||d_(6q)}dt 

to  to 

and  finally  be  appropriate  substitution  and  rearrangement. 
The  expression  for  variation  in  u  becomes: 

«u-  t|f««l*'+  /'{f  -  ffeC|§>>«>  «. 

"to   ^0 

Since  by  boundary  conditions,  the  curve  q=q(t)  always  passes 
through  two  fixed  points  t=t 0  and  t=t j ,  it  follows  that  the 
6q=0  at  these  points.   Hence 

af   * : 
Hr»-5q]    =  0        and  so 

dq    . 

t  o 

(5)  6u  =  /'{jg-jj-cj*))*,  dt  . 

t  o 
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The  necessary  condition  for  the  u  to  have  an  extreme  value 
for  variation  in  the  function  q=q(t)  as  obtained  by  varying 
the  parameters  a,& ,y , . . .    is  that  6u=0  for  all  values  of  6a, 
SB,    &y ,    •••,  i.e.  for  all  values  of  6q  from  Eq .  (5).   This 
is  only  true  if 

(6)     !■£  -  5r(|J)  "  0 
v  '  3q    dtv9q' 

This  result,  known  as  Euler's  equation,  determines  that 

function  q=q(t)  which  maximizes  or  minimizes  the  value  of  u. 

It  is  a  differential  equation  which  must  be  solved  to  give 

the  function  q=q(t)  sought.   Euler's  equation,  however,  is 

only  a  necessary  condition  for  extreme  values  of  u.   For 

distinguishing  between  maximum  and  minimum  values,  a  general 

criterion  is  not  readily  obtainable,  in  simple  practical 

case,  it  can  usually  be  judged  from  general  reasoning. 
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APPENDIX  C 


OECD  Countries  used  in  the  Analysis 


1. 

Canada 

2. 

United  States 

3. 

Japan 

4. 

Australia  *** 

5. 

Austria 

6. 

Belgium 

7. 

Denmark  * 

8. 

Finland  * 

9. 

France 

10. 

Germany 

11. 

Greece  * 

12. 

Iceland  ** 

13. 

Ireland  * 

14. 

Italy 

15. 

Netherlands 

16. 

Norway  * 

17. 

Portugal  * 

18. 

Spain 

19. 

Sweden  * 

20. 

Switzerland  * 

21. 

Turkey  * 

22. 

United  Kingdom 

*   Not  used  in  Natural  Gas 
Analysis . 

**  Used  only  in  Aviation  Fuel, 
Motor  Gasoline,  Kerosine 
and  Gas  Oil  Analysis. 

***   Used  only  in  Energy/Capita 
Analysis  except  for  the 
years  1966,  1967. 
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